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ABSTRACT 

We  report  spatially  resolved  spectroscopy  of  both  components  of  the  low-mass  pre-main-sequence  binary  GV  Tau. 
High-resolution  spectroscopy  in  the  K  and  L  bands  is  used  to  characterize  the  stellar  properties  of  the  binary  and  to 
explore  the  nature  of  the  circumstellar  environment.  We  find  that  the  southern  component,  GV  Tau  S,  is  a  radial 
velocity  variable,  possibly  as  a  result  of  an  unseen  low-mass  companion.  The  strong  warm  gaseous  HCN  absorption 
reported  previously  by  Gibb  and  coworkers  toward  GV  Tau  S  was  not  present  during  the  epoch  of  our  observations. 
Instead,  we  detect  warm  (~500  K)  molecular  absorption  with  similar  properties  toward  the  northern  infrared  com¬ 
panion,  GV  Tau  N.  At  the  epoch  of  our  observations,  the  absorbing  gas  toward  GV  Tau  N  was  approximately  at  the 
radial  velocity  of  the  GV  Tau  molecular  envelope,  but  it  was  redshifted  with  respect  to  the  star  by  ~  1 3  km  s  1 .  One 
interpretation  of  our  results  is  that  GV  Tau  N  is  also  a  binary  and  that  most  of  the  warm  molecular  absorption  arises  in 
a  circumbinary  disk  viewed  close  to  edge-on. 

Subject  headings:  circumstellar  matter  —  infrared:  stars  —  stars:  formation  —  stars:  fundamental  parameters  — 
stars:  individual  (GV  Tau)  —  stars:  late-type  —  stars:  low-mass,  brown  dwarfs  — 
stars:  pre-main-sequence  —  techniques:  radial  velocities  —  techniques:  spectroscopic 

Online  material:  color  figures 


1.  INTRODUCTION 

GV  Tau  (Haro  6-10,  Elias  7,  HBC  389,  IRAS  04263+2426)  is 
an  unusual  young  stellar  object  (YSO)  in  the  Taurus  molecular 
cloud.  While  the  spectral  energy  distribution  of  the  system  has 
the  strongly  rising  2-25  /im  continuum  that  is  characteristic  of 
Class  I  sources  (e.g.,  Furlan  et  al.  2007),  millimeter  studies  find 
comparatively  weak  dust  and  gas  emission  from  the  source  com¬ 
pared  to  other  low-mass  embedded  YSOs  (Hogerheijde  et  al. 
1998).  This,  and  the  poorly  defined  molecular  outflow  structure 
from  the  source,  suggest  that  the  system  lacks  a  significant  en¬ 
velope  component  (Hogerheijde  et  al.  1998)  and  may  be  a  more 
evolved  Class  I  system. 

GV  Tau  is  also  a  binary  (projected  separation  1.2";  Leinert  & 
Haas  1989).  It  is  one  of  a  small  population  of  low-mass  pre- 
main-sequence  binaries  in  which  one  component  is  optically  vis¬ 
ible,  while  the  other  is  optically  faint  and  radiates  primarily  at 
infrared  wavelengths.  The  optically  visible  component  in  the  GV 
Tau  system  is  the  southern  component  (GV  Tau  S),  a  T  Tauri  star. 
Its  stellar  properties  have  been  characterized  previously  based  on 
high-resolution  infrared  (Doppmann  et  al.  2005)  and  optical 
(White  &  Hillenbrand  2004)  spectroscopy.  GV  Tau  S  is  unusual 
in  that  its  radial  velocity  has  been  found  to  differ  significantly  (by 
more  than  3  a)  from  both  the  radial  velocity  of  the  surrounding 
molecular  cloud  (by  >7.5  km  s~')  and  the  radial  velocity  distri¬ 
bution  of  other  Taurus  sources  (Doppmann  et  al.  2005;  Covey  et  al. 
2006).  In  addition,  L-  and  M- band  absorption  by  gaseous  warm 
organic  molecules  ( HCN,  C2H2,  CO)  has  been  reported  toward 

1  Data  presented  herein  were  obtained  at  the  W.  M.  Keck  Observatory  from 
telescope  time  allocated  to  the  National  Aeronautics  and  Space  Administration 
through  the  agency’s  scientific  partnership  with  the  California  Institute  of  Tech¬ 
nology  and  the  University  of  California.  The  Observatory  was  made  possible  by 
the  generous  financial  support  of  the  W.  M.  Keck  Foundation. 

2  NOAO,  950  North  Cherry  Avenue,  Tucson,  AZ  85719;  gdoppmann@noao 
.edu,  najita@noao.edu. 

3  Naval  Research  Laboratory,  Code  7210,  Washington,  DC  20375;  carr@ 
nrl.navy.mil. 


GV  Tau  S  (Gibb  et  al.  2007),  making  it  one  of  the  few  pre-main- 
sequence  sources  in  which  such  absorption  features  have  been 
detected. 

Less  is  known  about  the  northern  component,  the  infrared  com¬ 
panion  of  the  binary  system  (GV  Tau  N ).  It  is  approximately  1 00  times 
fainter  than  G  V  Tau  S  at  optical  wavelengths  ( K.  Stapelfeldt  et  al. , 
in  preparation),  is  bright  at  infrared  wavelengths,  and  experiences 
significant  extinction,  as  indicated  by  the  strong  water  ice  and 
silicate  absorption  observed  toward  the  source  ( Whittet  et  al.  1988; 
Van  Cleve  et  al.  1994;  Leinert  et  al.  2001).  Extended  emission  (on  a 
10"  scale)  associated  with  the  GV  Tau  system  is  found  to  be  highly 
polarized  (6%)  at  infrared  wavelengths,  a  result  that  is  attributed  to 
scattering  by  a  flattened  envelope  (i.e.,  a  shell  with  an  aspect  ratio  of 
10: 1)  or  a  disk  viewed  at  high  inclination  (Menard  et  al.  1993). 
Since  strong  silicate  absorption  is  found  only  in  the  spectrum  of 
GV  Tau  N  and  not  GV  Tau  S,  this  is  taken  as  evidence  for  an 
edge-on  disk  associated  with  GV  Tau  N  (Van  Cleve  et  al.  1994). 

High  angular  resolution  observations  provide  limited  support 
for  this  interpretation.  When  studied  at  0.08"  resolution  at  K,  GV 
Tau  N  is  found  to  be  surrounded  by  a  nebula  whose  shape  may  be 
consistent  with  the  presence  of  an  edge-on  disk,  but  the  shape  of 
the  nebula  is  irregular  and  varies  with  time  ( Koresko  et  al.  1 999). 
The  irregular  structure  is  in  marked  contrast  to  the  symmetric 
lenslike  shape  of  scattered  light  distribution  seen  in  the  HK  Tau 
B  system,  a  system  that  is  also  believed  to  possess  a  nearly  edge-on 
disk  (Stapelfeldt  et  al.  1998;  Koresko  1998).  Adding  further  com¬ 
plexity,  significant  photometric  variability  is  found  for  both  com¬ 
ponents  at  1 .6-4.5  gm  on  timescales  as  short  as  a  month  (Leinert 
et  al.  2001),  with  the  3.1  gm  water  ice  band  absorption  strength 
also  varying  toward  both  components. 

Here  we  use  spatially  resolved  high-resolution  infrared  spec¬ 
troscopy  of  both  binary  components  to  characterize  the  stellar  and 
circumstellar  components  of  the  GV  Tau  system.  Our  K-  and 
L-band  spectroscopic  observations  are  described  in  §  2  and  the 
detected  spectral  features  are  described  in  §  3.  We  use  these  data 
to  characterize  the  stellar  and  circumstellar  properties  of  the 
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system  (§  4).  The  results  are  discussed  in  §  5  and  our  conclusions 
are  presented  in  §  6. 

2.  OBSERVATIONS  AND  DATA  REDUCTION 
2.1.  Spectroscopic  Observations 

K-  and  Z-band  infrared  spectra  were  obtained  on  2007  January  3 
using  the  cryogenic  echelle  facility  spectrometer,  NIRSPEC 
(McLean  et  al.  1998),  on  the  10  m  Keck  II  telescope  atop  Mauna 
Kea,  Hawaii.  Spectra  were  acquired  through  the  0.432"  (3  pixel) 
wide  slit,  providing  spectroscopic  resolution  R  =  )J6).  =  24,000 
(12.5  km  s_1).  The  echelle  and  cross  disperser  angles  were  ori¬ 
ented  to  obtain  Z-band  spectra  of  transitions  of  HCN  and  C2H2 
(e.g.,  3.04-3.09  pm,  order  25).  We  used  the  KL  blocking  filter  to 
image  the  Z-band  orders  onto  the  instrument’s  1024  x  1024  InSb 
detector.  Immediately  following  the  Z-band  observations,  the 
echelle  and  cross  disperser  gratings  were  reoriented  to  allow  key 
A'- hand  features  (i.e.,  Mg  1,  Al  1,  Na  1, ;/  =  2-0  CO,  and  Br7)  to 
fall  onto  the  detector.  In  this  configuration,  orders  32-38  (non¬ 
contiguous)  were  imaged  through  the  NIRSPEC-7  blocking  fil¬ 
ter  onto  the  detector  array. 

Both  binary  components  of  GV  Tau  were  acquired  in  the  slit 
simultaneously,  at  a  position  angle  on  the  sky  of  356°  (Leinert  & 
Haas  1989).  Initially,  GV  Tau  N  (the  brighter  component  in  our 
observations  of  X  >  2.09  pm)  was  centered  in  the  slit.  During 
the  Z-  and  K-band  observations  the  slit  was  held  physically  sta¬ 
tionary  to  avoid  slight  grating  angle  shifts  caused  by  vibrations 
from  the  slit  rotator  motor.  During  the  time  we  integrated  on  the 
binary  pair  (<25  minutes  total),  the  measured  flux  ratio  of  the  com¬ 
ponents  in  the  slit  did  not  change  systematically.  This  is  consistent 
with  the  negligible  expected  motion  of  the  binary  pair  within  the  slit 
as  the  nonequatorially  mounted  telescope  tracked.  Despite  some 
wind  shake  in  the  direction  across  the  slit  in  the  east-west  direction, 
the  individual  components  (separated  by  1 .2")  were  well  resolved 
in  median  seeing  (^0.6"),  allowing  separate  spectral  extractions 
of  the  southern  and  northern  components. 

The  data  were  acquired  in  pairs  of  exposures  of  40  and  60  s  in 
the  Z  and  A  bands,  respectively.  The  telescope  was  nodded  ±6" 
from  the  slit  center  in  an  ABBA  sequence  along  the  24"  long  slit, 
so  that  each  frame  pair  would  contain  object  and  sky  in  both  nod 
positions.  To  correct  for  telluric  absorption,  we  also  obtained  spec¬ 
tra  of  an  early-type  star  located  close  to  GV  Tau  in  the  sky  (HR 
1412),  before  and  after  the  GV  Tau  observations  (Z-  and  A'-band 
settings,  respectively).  Spectra  of  the  internal  NIRSPEC  contin¬ 
uum  lamp  were  taken  for  flat  fields  at  the  K-  and  Z-band  settings. 
Exposures  of  the  argon,  neon,  krypton,  and  xenon  arc  lamps  pro¬ 
vided  wavelength  calibration  for  all  the  K-band  orders,  except  in 
order  33,  where  we  used  the  telluric  lines  for  calibration  due  to  the 
scarcity  of  arc  lamp  lines  in  this  order.  We  also  obtained  short  ex¬ 
posure  (10x0.1  s)  K-  and  Z-band  images  of  the  GV  Tau  binary 
using  the  slit  viewing  camera  detector  (SCAM )  in  order  to  mea¬ 
sure  the  relative  brightness  of  the  two  components. 

2.2.  Data  Reduction 

The  spectra  were  reduced  using  standard  IRAF  packages  ( Massey 
et  al.  1992;  Massey  1997).  Sky-subtracted  beam  pairs  were  divided 
by  a  normalized  flat  field.  Individual  echelle  orders  were  parsed 
from  the  multiorder,  co-added  images.  Bad  pixels  (i.e.,  hot  or 
low-responsivity  pixels)  in  each  order  were  identified  by  inspec¬ 
tion,  and  removed  by  interpolation  using  f  ixpix  in  IRAF.  To 
better  remove  telluric  emission  features  which  are  more  severe  in 
the  Z  band,  we  rectified  the  slit  in  our  Z-band  data  using  a  third- 
order  transform  solution  derived  from  the  brightest  emission  sky 
lines  (typically  10-12  per  order)  traced  along  the  slit  length  at 


each  beam  position  before  sky  subtraction.  Transfomied  and 
cleaned  orders  were  extracted  using  the  apall  task.  Extractions 
for  each  binary  component  of  GV  Tau  were  based  on  a  signal 
profile  down  to  50%  (FWHM)  on  both  sides  of  the  profile  peak, 
to  minimize  the  contaminating  flux  from  the  nearby  companion. 
For  the  Z-band  spectra,  residual  background  sky  was  subtracted 
in  the  extraction  using  selected  regions  of  the  profile  cut  that 
were  well  outside  of  the  profile  of  the  spatially  double-peaked 
binary. 

Wavelength  calibration  in  the  Z  band  was  achieved  using 
selected  telluric  absorption  features  present  in  the  spectra  of  the 
telluric  standard  star  ( HR  1412).  Stronger  unblended  telluric  lines 
were  selected  over  weaker  ones,  while  avoiding  very  saturated 
lines  whose  depths  went  to  zero.  Rest  wavelengths  for  the  absorp¬ 
tion  features  were  obtained  from  the  H1TRAN  database  (  Rothman 
et  al.  1998).  Most  or  all  of  these  same  lines  were  used  to  derive  a 
wavelength  solution  for  the  target,  with  the  exception  of  those  lines 
that  might  fall  very  close  in  wavelength  to  emission  or  absorption 
lines  present  in  the  target  spectrum. 

Telluric  features  present  in  the  spectra  of  GV  Tau  at  each  nod 
position  in  the  slit  were  removed  by  dividing  by  the  spectrum  of 
the  telluric  standard  (HR  1412)  observed  at  the  same  nod  po¬ 
sition.  Several  weak  stellar  absorption  lines  were  present  in  the 
K-band  orders  of  the  standard  due  to  its  relatively  late  spectral 
type  (A7  III).  These  atomic  lines  were  modeled  and  removed  by 
dividing  the  standard  spectrum  by  the  best-fit  stellar  synthetic 
spectrum  before  division  into  the  GV  Tau  spectrum.  The  stellar 
synthetic  spectra  were  generated  with  the  program  MOOG  (Sneden 
1973)  using  the  NextGen  model  atmospheres  (Hauschildt  et  al. 
1999).  The  initial  line  list  was  taken  from  Kurucz  (1993)  and  in¬ 
dividual  line  parameters  were  adjusted  to  fit  the  observed  disk- 
center  solar  spectrum  of  Livingston  &  Wallace  (1 99 1 ).  For  HR  1412, 
we  used  a  model  with  Zeff  =  7600  K  and  log  g  =  3.5.  Small  ad¬ 
justments  were  made  in  the  elemental  abundances  to  give  the  best 
empirical  fit  to  HR  1412.  We  measured  t>sin/  =  75  km  s-1,  in 
agreement  with  published  values  (Royer  et  al.  2002). 

Br7  absorption  at  2.166  pm  was  also  present  in  the  telluric 
standard.  We  therefore  used  the  synthetic  telluric  spectral  mod¬ 
eling  program,  AT  (E.  Grossman  1989,  private  communication), 
to  model  the  telluric  features  that  we  observed  in  this  region  of 
the  HR  1412  spectrum.  We  fit  the  depths  of  the  observed  telluric 
lines  and  then  used  the  best-fit  synthetic  spectrum  to  divide  out 
the  telluric  lines  in  the  GV  Tau  spectra  in  this  order  (order  35, 
Fig.  lc).  In  order  to  view  the  true  continuum  shape  within  each 
echelle  order,  we  then  multiplied  the  resulting  GV  Tau  spectra  by 
the  spectral  slope  of  a  7600  K  blackbody  (Figs,  la-le). 

Simple  aperture  photometry  of  the  SCAM  images  showed  that 
GV  Tau  N  was  35%-40%  brighter  than  GV  Tau  S  in  the  K+open, 
N1RSPEC7,  and  KL  filters.  The  relative  flux  of  the  GV  Tau  N  and 
S  K-band  spectra  in  orders  33  and  34  (Figs,  la  andlb),  which  are 
located  in  the  middle  of  the  K-band  filters  above,  are  consistent 
with  the  aperture  photometry. 

3.  DETECTED  SPECTRAL  FEATURES 

Emission  and  absorption  features  in  GV  Tau  N  and  S  were 
detected  in  the  K-band  spectra  of  each  object  (Figs,  la-ld).  In 
order  33  (2.269-2.304  pm;  Fig.  la),  v  =  2-0  CO  overtone  emis¬ 
sion  and  absorption  are  present  in  both  GV  Tau  N  and  S.  The 
emission  component  dominates  the  spectrum  in  GV  Tau  N  in  this 
order,  while  the  absorption  dominates  in  GV  Tau  S.  Lines  due  to 
absorption  of  neutral  atomic  species  that  are  present  in  GV  Tau  S 
in  this  order  (e.g.,  Fe  1  at  2.2747  pm  and  Mg  1  at  2.2814  pm)  are 
absent  in  GV  Tau  N.  The  A'-band  absorption  features  in  both  com¬ 
ponents  are  presumed  to  originate  from  the  stellar  photosphere. 
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Fig.  1  a 

Fig.  1 . —  Simultaneous  spectra  of  GV  Tau  N  (top)  and  GV  Tau  S  (bottom)  in  the  band  in  (a)  order  33,  (b)  order  34,  (c)  order  35,  (d)  order  36,  and  (e)  in  the Z-band  or¬ 
der  25.  The  relative  flux  ratio  ( Aouth//North)  of  the  two  components  is  (a)  ~0.7  (order  33),  (b)  0.74  (order  34),  (c)  0.85  (order  35),  (d)  ~  1 .0  (order  36),  and  (e)  0.45  (order  25). 
The  A'-band  orders  show  emission  features  in  both  objects,  while  absorption  lines  dominate  the  spectra  in  GV  Tau  S.  (e)  The  A-band  spectrum  of  GV  Tau  N  shows  strong 
P-branch  absorption  lines  of  HCN  (gray  tick  marks,  labeled)  and  several  weaker  C2H2  lines  (black  tick  marks),  while  these  features  are  absent  in  GV  Tau  S.  Spectral  regions 
with  poor  telluric  transmission  have  been  excised  from  this  panel.  [See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this  figure .] 


Emission  features  of  Na  i  and  EE  dominate  the  structure  of  the 
GV  Tau  N  spectrum  in  order  34  (2.203-2.236  /im;  Fig.  lb).  In¬ 
terestingly,  the  broad  emission  lines  near  2.2065  and  2.2090  /im 
in  GV  Tau  N  appear  to  show  absorption  components  possibly  of 
a  photospheric  origin.  The  £(0)1—0  H2  emission  at  2.2235  /im  is 
detected  in  both  components,  and  is  spatially  extended  (see  §  4.3). 
Other  than  the  H2  emission,  the  spectrum  of  GV  Tau  S  in  this  or¬ 
der  displays  neutral  atomic  absorption  lines  of  photospheric  ori¬ 
gin  (Fig.  lb). 

The  spectra  in  order  35  (2.140-2.173  /im,  Fig.  lc)  are  char¬ 
acterized  by  strong  H  i  Br7  emission  in  both  GV  Tau  N  and  S. 
The  continuum  level  near  the  Brackett  7  line  is  ~20%  higher  in 
GV  Tau  N,  while  the  equivalent  width  of  the  emission  line  is  3  times 
greater  in  GV  Tau  S. 

In  order  36  (2.081-2.113  /im;  Fig.  Id),  the  spectrum  of  GV 
Tau  N  appears  featureless.  In  contrast,  GV  Tau  S  exhibits  pho¬ 
tospheric  absorption  lines  (e.g..  Mg  1  at  2.1065  /im  and  A1 1  at 
2.1099  /im),  which  are  useful  as  diagnostics  of  the  stellar  effec¬ 
tive  temperature  and  surface  gravity. 


With  ourZ-band  spectra  in  order  25  (3.0450-3.0865  /an;  Fig.  le), 
we  detect  strong  HCN  absorption  (~10%  deep)  and  weaker 
C2H2  absorption  in  GV  Tau  N.  In  contrast,  no  molecular  absorp¬ 
tion  is  detected  in  the  GV  Tau  S  spectrum.  In  §  4.4  we  charac¬ 
terize  the  molecular  absorption  in  GV  Tau  N. 

4.  RESULTS 

4.1.  Stellar  and  Circumstellar  Properties 
4.1.1.  South  Component 

All  four  Zl-band  orders  of  the  GV  Tau  S  spectrum  display 
stellar  photospheric  absorption  by  neutral  atomic  features  (e.g., 
Na  1,  Si  1,  Sc  1,  Mg  1,  A1 1)  and  v  =  2-0  CO.  We  used  synthetic 
stellar  spectral  models  generated  by  the  synthesis  program 
MOOG  (Sneden  1973)  and  using  NextGen  stellar  atmosphere 
model  structures  (Hauschildt  et  al.  1999)  to  constrain  the  stellar 
properties  of  GV  Tau  S.  Our  modeling  focused  on  three  spectral 
subregions:  The  “Na”  region  (2.204-2.210  /tm),  the  “Mg/Al” 
region  (2. 104-2.  Ill  /im),  and  the  “CO”  region  (2.292-2.300  /im). 
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Fig.  1  b 


These  regions  of  the  A'-band  contain  the  strongest  absorption  fea¬ 
tures  with  which  to  measure  the  stellar  properties  of  low-mass 
YSOs  (Doppmann  et  al.  2005).  The  absorption  lines  of  Na  i  and 
Mg  i  are  particularly  gravity  and  temperature  sensitive,  but  in  the 
opposite  sense  from  each  other.  For  example,  at  cool  effective 
temperatures  (3200-4500  K)  and  subdwarf  surface  gravities 
(3.5  <  log  g  <  4.5)  Na  and  Mg  lines  both  grow  stronger  as  log  g 
increases,  but  an  increase  in  causes  Mg  lines  to  grow  while 
the  Na  lines  weaken.  A  simultaneous  fit  of  synthesis  models  to  the 
Na  and  Mg/Al  regions  (Fig.  2)  constrains  these  key  stellar  prop¬ 
erties  in  reducing  or  breaking  the  temperature-gravity  degeneracy 
displayed  in  a  single  absorbing  species  (Doppmann  et  al.  2005). 

The  best  model  fit  to  these  wavelength  subregions  yields 
Teff  =  3800  ±  100  K,  logyy  =  4.0  ±  0.17,  tisin;  (rotation)  = 
24  ±  3  km  s-1,  and  A-band  veiling  rK  =  2.5  ±  0.2  (Fig.  2). 
From  the  Teg  and  logyy  derived  above,  we  estimate  a  stellar 
luminosity,  mass,  and  radius  of  C*  =  0.3  £0,Af»  =  0.5  M@,  and 
R*  =  1.2  R  .  using  the  pre-main-sequence  model  tracks  of  Siess 
et  al.  (2000),  where  a  3  Myr  isochrone  was  consistent  with  our 
values  of  temperature  and  gravity.  The  derived  luminosity  and 
effective  temperature  imply  an  estimated  K- band  extinction  to¬ 
ward  GV  Tau  S  of  Ak  =  0.7  based  on  the  most  recent  photo- 
metery  from  the  literature  (2000  March  7;  Leinert  et  al.  2001), 


which  reported  a  value  (m^  =  8.61)  that  was  within  1  cr  of  the 
average  in  the  study  over  the  time  period  from  1988  September 
to  2000  March  (mjc  =  8.1  ±  0.7).  Our  derived  K- band  extinc¬ 
tion  is  consistent  with  the  average  optical  extinction  of  A  v  =  5.6 
estimated  by  Koresko  et  al.  (1997). 

Our  derived  Teg,  v  sin  i,  and  j>Lsr  values  are  consistent  with 
the  high-resolution  optical  study  of  White  &  Hillenbrand 
(2004)  given  the  low  signal-to-noise  ratio  (S/N  <  10)  in  their 
1999  December  6  observation  of  GV  Tau  S.  The  moderately 
high  veiling  we  derive  (r^  =  2.5)  is  more  characteristic  of  Class  1 
or  flat-spectrum  YSOs  than  of  T  Tauri  stars,  although  still  within 
the  range  of  veiling  values  measured  in  a  sample  of  Ophiuchus 
Class  II  sources  (Doppmann  et  al.  2003). 

These  results  agree  well  with  the  veiling  measured  from  the  pre¬ 
vious  observation  of  GV  Tau  S  from  2001  November  6  (Doppmann 
et  al.  2005).  Doppmann  et  al.  (2005)  quote  a  veiling  value  of 
i’k  =1.8,  which  includes  a  correction  for  a  systematic  effect  seen 
in  the  best-fit  synthesis  models  to  observations  of  MK  standards 
(see  eq.  [1]  in  §  3.6  of  Doppmann  et  al.  2005).  The  measured  veil¬ 
ing  without  the  correction  for  the  systematic  effect  was  i~k  =  2.4, 
which  is  consistent  with  our  measured  value. 

The  constraint  on  effective  temperature  and  surface  gravity  is  con¬ 
sistent  with  the  results  from  Doppmann  et  al.  (2005).  Doppmann 
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et  al.  (2005)  reported  a  somewhat  warmer  effective  temperature 
(Teff  =  4500  K)  and  higher  surface  gravity  (log  g  =  4.4),  which 
was  based  on  fits  to  synthetic  spectra  of  a  smaller  number  of 
photospheric  absorption  features  (excluding  Sc  i  lines  at  2.2058 
and  2.207 1  gm)  at  a  similar  signal-to-noise  ratio  and  slightly  lower 
spectral  resolution  ( R  =  18,000).  A  comparably  good  fit  to  the 
2001  data  can  be  obtained  with  =  4000  K  and  log  g  =  3.9 
that  includes  the  Sc  i  lines  above  as  was  done  in  this  study. 

The  radial  velocity  of  GV  Tau  S  was  measured  by  shifting  a 
synthetic  template  spectrum  to  independently  match  each  of  12 
selected  absorption  lines  that  were  present  in  orders  33  and  34 
(see  Table  1).  The  error  was  determined  from  the  quadrature  sum 
of  the  systematic  error,  dominated  by  the  uncertainty  in  the  wave¬ 
length  calibration  of  the  data,  and  the  statistical  error,  measured 
from  the  variation  in  the  radial  velocity  values  among  the  indi¬ 
vidual  absorption  lines.  Our  measured  im.sr  of  GV  Tau  S  differs 
from  previous  measurements  in  the  literature  for  this  source.  In 
2001,  the  radial  velocity  of  GV  Tau  S  was  measured  to  be  j’i.sr  = 
—6.2  ±1.5  km  s_1  (Covey  et  al.  2006),  blueshifted  by  ~9- 
1 5  km  s~ 1  relative  to  the  observations  of  this  study  and  the  value  of 
flsr  =  ±3.1  ±  3.8  km  s_I  measured  by  White  &  Hillenbrand 
(2004).  As  a  result,  this  source  was  flagged  as  a  radial  velocity  out¬ 


lier,  since  its  radial  velocity  was  >3  er  from  the  mean  of  all  the 
Taurus  sources  in  the  survey  (Covey  et  al.  2006).  Our  2007  obser¬ 
vations  now  place  the  radial  velocity  of  GV  Tau  S  (i>lsr  =  9.4  ± 
1.7  km  s~')  close  to  the  mean  olsr  (4.6  ±1.8  km  s~')  of  the 
sources  in  the  2001  Taurus  survey  (Doppmann  et  al.  2005;  Covey 
et  al.  2006).  We  discuss  the  radial  velocity  variations  further  in  §  5 . 1 . 

As  described  in  §  3,  the  CO  v  =  2-0  bandhead  of  GV  Tau  S 
has  both  emission  and  absorption  components.  Using  the  best-fit 
stellar  parameters  obtained  from  the  Na  and  Mg/Al  regions  of 
the  GV  Tau  S  spectrum  as  a  constraint  on  the  stellar  photospheric 
properties,  we  fit  the  CO  feature  with  a  composite  model  of 
circumstellar  disk  emission  and  stellar  photospheric  absorption 
(Fig.  3).  The  emission  component  is  modeled  as  arising  from  a 
differentially  rotating  disk  that  is  in  chemical  equilibrium.  The 
emission  arises  between  an  inner  radius  Rj,  at  which  the  pro¬ 
jected  disk  rotational  velocity  has  a  specified  value  of  v  sin  i,  and 
an  outer  radius  R0.  The  disk  temperature  and  column  density  are 
modeled  as  simple  radial  gradients  that  decrease  as  power  laws 
(Carr  1993;  Carr  et  al.  2004).  The  star  and  disk  are  placed  at  the 
same  radial  velocity  and  the  composite  model  spectrum  is  smoothed 
to  the  13  km  s  1  velocity  resolution  of  our  NIRSPEC  data.  The 
model  fit  parameters  are  given  in  the  caption  to  Figure  3. 
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4.1.2.  North  Component 

The  spectrum  of  GV  Tail  N  shows  emission  and  absorption 
components  in  both  the  CO  and  Na  regions  (Figs.  4  and  5)  of  the 
K  band.  We  attribute  the  emission  to  a  circumstellar  disk;  the 
absorption  features  are  presumed  to  arise  from  the  stellar  photo¬ 
sphere,  implying  a  late  spectral  type  for  the  star.  We  assume  a  stel¬ 
lar  mass  of  M*  =  0.8  Me  and  an  age  of  r  =  3  Myr  for  this  object, 
consistent  with  the  typical  values  found  for  low-mass  YSOs  in  Taurus 
( Briceno  et  al.  2002).  These  assumptions  imply  an  effective  tem¬ 
perature  of  Teff  =  4 1 00  K  and  a  surface  gravity  of  log  g  ^  4 . 0  from 
pre-main-sequence  evolutionary  model  tracks  (Siess  et  al.  2000). 

Using  these  assumed  stellar  properties,  we  modeled  the  CO 
feature  observed  in  GV  Tau  N  as  a  combination  of  circumstellar 
disk  emission  and  stellar  photospheric  absorption.  As  in  the  case 
of  GV  Tau  S,  the  emission  component  is  fit  with  a  simple  model 
of  CO  emission  from  a  differentially  rotating  disk.  The  details  of 
the  model  fit  are  given  in  the  caption  to  Figure  4.  We  find  a  good  fit 
to  the  observed  CO  overtone  emission  ( Fig.  4)  with  the  additional 
assumptions  of  slow  stellar  rotation  (i)sin/  =  15  km  s_1)  and 
heavy  veiling  (tp  =  12).  Such  a  substantial  infrared  veiling  greatly 
exceeds  what  has  been  measured  from  near-IR  photospheric  ab¬ 
sorption  lines  in  past  studies  (e.g.,  rK  <  4.5;  Luhrnan  et  al.  1998; 


Luhman  &  Rieke  1999;  Doppmann  et  al.  2003, 2005).  It  also  results 
in  large  uncertainties  in  the  inferred  W-band  extinction. 

Broad  Na  emission  has  been  observed  in  several  Class  I  and 
flat-spectrum  protostars,  always  accompanied  by  v  =  2-0  CO 
emission  ( Doppmann  et  al.  2005).  The  absorption  lines  in  the  Na 
region,  apparent  within  the  broader  emission  features,  can  be  fit 
with  the  same  stellar  photospheric  model  and  somewhat  larger 
veiling  (rK  =  15;  Fig.  5)  than  was  used  in  fitting  the  CO  region. 
If  we  subtract  the  stellar  photospheric  model  from  the  observed 
spectrum,  the  resulting  spectrum  is  similar  to  the  Na  emis¬ 
sion  seen  in  other  T  Tauri  stars.  The  absence  of  detected  spec¬ 
tral  features  in  the  Mg/Al  region  of  GV  Tau  N  in  our  data 
(signal-to-noise  ratio  ~  170)  is  also  consistent  with  the  model 
assumptions. 

The  radial  velocities  of  the  CO  and  Na  absorption  compo¬ 
nents  in  GV  Tau  N  are  in  agreement  with  one  another  (i>lsr  = 
—4.5  km  s-1;  see  Table  2),  and  with  the  radial  velocity  of  the 
emission  component  from  which  they  arise.  Given  the  assumed 
mass  M*  =  0.8  Mg,  we  estimate  the  stellar  luminosity,  radius, 
and  A'- band  extinction  (i.e.,  L,  =  0.6  LQ,  R,  1.4  Re,  and 
Ak  =  2.6)  from  the  photospheric  absorption  we  find  in  our  spec¬ 
tra,  using  Siess  et  al.  (2000)  evolutionary  model  tracks  and  Leinert 
et  al.  (2001)  photometry  from  2000  March  7  (mx  =  8.66). 
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The  unknown  apparent  magnitude  of  GV  Tau  N  on  the  date  of 
our  observations  introduces  an  uncertainty  in  the  derived  K-band 
extinction  and  veiling.  If  the  apparent  magnitude  was  the  his¬ 
torical  average  of  mx  =  9.8  (Leinert  et  al.  2001),  instead  of  the 
value  reported  on  2000  March  7  of  mx  =  8.66  (Leinert  et  al. 
2001),  then  the  extinction  (Ax)  would  increase  by  1.2.  If  we 
had  assumed  a  higher  (Teff  =  4300  K)  or  lower  (Teff  =  3900  K) 
temperature  than  the  value  we  use  here  (Teff  =  4100  K),  which 
was  based  on  the  IMF  peak  (Briceno  et  al.  2002),  this  would 
imply  a  veiling  (rx)  of  1 1  or  13,  and  an  extinction  (A Ax)  of  +0.5 
or  —0.6  relative  to  the  derived  value  presented  above  (AK  =  2.6). 

4.2.  Brackett  y  Emission 

We  detect  H  1  Br7  emission  in  both  components  of  GV  Tau 
(Fig.  6).  The  broad  emission  profiles  (i.e.,  175  and  170  km  s_1 
FWFLM,  GV  Tau  S  and  N,  respectively)  are  consistent  with  the 
Br7  width  of  other  active  T  Tauri  stars  (Najita  et  al.  1996;  Folha 
&  Emerson  2001).  In  GV  Tau  S,  where  the  equivalent  width  is  3 
times  greater  than  in  GV  Tau  N  (EWsouth  =  —  3.9  A),  the  cen¬ 
troid  velocity  is  very  close  to  the  stellar  velocity,  as  measured  by 
the  stellar  photospheric  absorption  lines  (see  §  4.1.1).  However, 
the  centroid  velocity  of  the  Br7  line  in  GV  Tau  N  is  blueshifted 
by  12.5  km  s_1  relative  to  the  CO  absorption  component  (§  4. 1 .2). 


Blueshifted  Bry  centroids,  which  are  common  among  T  Tauri 
stars,  are  consistent  with  an  origin  for  the  emission  in  gas  infalling 
in  a  stellar  magnetosphere  (Najita  et  al.  1996;  Muzerolle  et  al. 
1998b;  Folha  &  Emerson  2001). 

We  can  use  the  Bry  line  strength  to  place  a  rough  constraint  on 
the  contribution  of  stellar  accretion  in  GV  Tau  N  and  S  to  the 
bolometric  luminosity  of  the  system.  We  first  convert  the  emis¬ 
sion  equivalent  width  to  a  line  luminosity  using  the  2000  March  7 
K-band  photometry  from  Leinert  et  al.  (2001)  and  the  estimated 
K-band  extinction  for  each  component  (§§4.1.2  and  4.1.1).  The 
estimated  Bry  luminosities  for  GV  Tau  N  and  S  are  9  x  1 0  5  and 
5  x  10-5  Lq  ,  respectively.  Using  the  empirical  relation  given  by 
Muzerolle  et  al.  (1998b)  in  their  Figure  4,  we  infer  accretion 
luminosities  of  0. 1  and  0.2  L@  for  GV  Tau  N  and  S,  respectively. 
Thus,  we  estimate  the  total  luminosity  (stellar+hot  accretion)  of 
the  GV  Tau  system  to  be  1.2  ZQ,  below  what  has  been  estimated 
for  its  bolometric  system  luminosity  (Zboi  =  7-9  Z0;  Kenyon  & 
Hartmann  1995;  Furlan  et  al.  2007). 

4.3.  H2  Emission 

The  5(0)  1-0  H2  emission  (2.22329  /.tm)  that  we  detect  from 
GV  Tau  is  marginally  resolved  in  our  spectra  with  a  velocity  width 
of  14  and  17  km  s_1  for  GV  Tau  S  and  N,  respectively.  Spatially, 
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Fig.  2. — Two  regions  of  the  K  band  of  GV  Tau  S  (histogram)  that  show  stellar  photospheric  absorption  lines.  A  spectral  synthesis  model  (dotted  curve)  with 
reff  =  3800  K,  log  g  =  4.0,  tisin i  =  24  km  s-1  (includes  instrumental  broadening),  along  with  a  veiling  of  vk  —  2.5  fits  the  spectrum.  The  model  fit  in  the  top  panel  is  at  a 
radial  velocity  of  vlsr  =  9.1  km  s-1  (see  Table  1).  [See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this  figure.] 


the  bulk  of  the  emission  is  coincident  with  GV  Tau  N  and  S 
(within  the  0.5"  extraction  aperture)  but  extends  along  the  slit 
to  the  south  ( P.A.  =  176°)  ending  at  a  bright  knot  6.5"  south  of 
GV  Tau  S  (Fig.  7). 

HST  imaging  of  the  GV  Tau  system  shows  resolved  nebular 
emission  extending  away  in  an  arc  to  the  east,  south,  and  southwest 
of  GV  Tau  S  while  GV  Tau  N  is  comparatively  free  of  nebulosity 
(K.  Stapelfeldt  2007,  private  communication).  The  location  of 
the  emission  knot  is  spatially  coincident  with  a  wisp  of  faint  op¬ 
tical  nebulosity  in  the  HST  image.  No  near-IR  line  emission  is  seen 
immediately  to  the  north  of  GV  Tau,  which  is  coincident  with  a 
lack  of  optical  nebulosity  in  this  region.  The  radial  velocity  of  the 
HE  emission  is  constant  along  the  slit  (plsr  =  8.5  km  s_1)  and 
agrees  with  the  radial  velocity  of  the  HCN  absorption  (§  4.4)  and 
molecular  envelope  for  this  system  (Hogerheijde  et  al.  1998). 

In  order  37,  we  detect  5(2)1 -0  H2  emission  in  both  stellar 
components.  The  emission  is  extended  southward  along  the  slit 
in  the  same  way  as  the  5(0)  emission.  The  emission  is  2-3  times 
brighter  than  the  5(0)  line  emission.  As  with  the  5(0)  line,  the 
5(2)  equivalent  width  is  ^2  times  stronger  in  GV  Tau  S  compared 
to  GV  Tau  N.  In  both  components,  the  equivalent  widths  and  ve¬ 


locity  widths  are  greater  than  their  counterparts  in  the  5(0)  lines. 
The  observed  HE  Fine  ratios  [i.e.,  1  — 05(2)/l  —05(0)]  inGVTauN 
and  S  are  consistent  with  shock  excitation,  similar  to  most  of  the 
classical  T  Tauri  stars  observed  by  Beck  et  al.  (2007). 

4.4.  Molecular  Absorption 

In  order  to  characterize  the  molecular  absorption  detected  in 
the  /.-band  spectrum  of  GV  Tau  N,  we  used  a  simple  model  of 
absorption  by  a  slab  at  a  single  temperature  and  column  den¬ 
sity.  We  adopted  the  HCN  line  list  from  the  HITRAN  database 
(Rothman  et  al.  1998).  The  data  are  well  fit  with  a  slab  temper¬ 
ature  T  =  550  K  and  column  density  Vhcn  =  1 .5  x  1017  cm  2, 
with  a  microturbulent  line  broadening  of  i>turb  =  3  km  s_1  and 
radial  velocity  «LsR  =  8.7  ±1.0  km  s_1  (Fig.  8).  The  radial  ve¬ 
locity  of  the  molecular  absorption  was  measured  using  selected 
P-branch  lines  ( Table  3)  in  regions  of  the  spectrum  that  had  good 
telluric  transmission  and  that  were  unblended  with  other  species 
(such  as  C2H2).  We  constructed  an  average  HCN  absorption 
profile  from  four  selected  HCN  absorption  transitions  ( PI  1 ,  P 1 4, 
PI 5,  and  PI 6)  to  which  we  fit  a  Gaussian  profile  (Fig.  9).  The  3  cr 
error  of  1.0  km  s_1  on  the  radial  velocity  was  estimated  by 


306 


DOPPMANN,  NAJITA,  &  CARR 


Vol.  685 


TABLE  1 

GV  Tau  South:  Radial  Velocity  Measurements 


^LSR 

Order/ Wavelength  Lines(s)  Absorption  Emission  (km  s-1)  Notes 


Order  37,  2.03  gm .  £(2)  1-0  H2  No  Yes  10.3  Line  centroid 

Order  36,  2.11  /im .  Mg  1,  A1  1  Yes  No  ...  Radial  velocity  poorly  constrained  in  this  region 

Order  35,  2.16 /im .  Br7  No  Yes  3.7  Line  centroid 

Order  34,  2.21  /im .  Na  1,  Ti  1,  Fe  1  Yes  No  9.1  ±  1.3  Mean  value  of  7  lines 

Order  34,  2.22  /im .  £(0)1-0  H2  No  Yes  8.5  Emission  extends  south  of  the  star 

Order  33,  2.29  /im .  Fe  1,  Ca  1,  Ti  1  Yes  No  10.0  ±  1.3  Mean  value  of  4  lines 

Order  33,  2.29  /im .  is  =  2-0  CO  Yes  Yes  . . .  Radial  velocity  poorly  constrained  in  this  region 

Order  25,  3.07  /im .  HCN,  C2H2  No  No 


comparing  the  Gaussian  fit  of  the  average  absorption  profile  to 
each  of  the  P-branch  line  fits  separately. 

5.  DISCUSSION 

To  summarize,  both  GV  Tau  N  and  S  appear  to  possess  slowly 
rotating  (t>sin/  =  15-24  km  s_1)  late  type  (^K7-M2)  photo¬ 
spheres.  While  the  stellar  properties  are  consistent  with  an  age  of 
~3  Myr  based  on  current  pre-main-sequence  stellar  evolutionary 
tracks,  the  age  estimate  is  uncertain  because  ages  inferred  from 
the  tracks  have  yet  to  be  observationally  verified.  We  also  find 
that  GV  Tau  S  is  a  radial  velocity  variable.  For  GV  Tau  N,  we 
find  that  the  system  shows  strong  CO  overtone  emission,  strong 
Al-band  veiling  (rK  =  12-15),  and  a  blueshifted  Br7  centroid. 
These  properties  are  consistent  with  a  young  star  undergoing  active 
accretion  in  a  disk  and  stellar  magnetosphere.  The  stellar  radial 
velocity  found  for  GV  Tau  N  (t>LSR  =  —4.5  ±  4.0  km  s~')  dif¬ 
fers  from  the  radial  velocity  of  the  GV  Tau  molecular  envelope 
(r’LSR  =  7.0  ±  0.5  km  s-1;  Hogerheijde  et  al.  1998)  and  the 
radial  velocity  of  the  warm  (^550  K)  HCN  absorption  detected 
intheZ,band(i>LSR  =  8.7  ±  1.0  km  s~ 1 ).  We  discuss  the  possible 
implications  of  these  results  in  the  rest  of  this  section. 


Topocentric  Wavelength  (pm) 

Fig.  3. — CO  overtone  spectral  region  of  GV  Tau  S  {black  histogram)  showing 
both  emission  and  absorption  components.  We  fit  the  observed  spectrum  using  a 
model  of  disk  emission  and  stellar  photospheric  absorption  components  at  the 
same  radial  velocity.  Our  disk  emission  model  {dashed  line)  has  the  proper¬ 
ties  vsini  (inner  edge)  =  88  km  s-1,  R0/Ri  =  7,  T  —  3500 (r/R;)-0'5  K,  E  = 
0.5(r/R;)-0'5  g  cm-2,  where  r  is  the  disk  radius.  The  stellar  model  {dotted  line)  has 
a  temperature  7^  =  3800  K,  gravity  log  g  =  4.0,  rotation  vsin  i  —  24  km  s”1, 
and  veiling  vk  —  2.6.  The  radial  velocity  of  the  model  fit  (ulsr  —  9.4  km  s-1)  is 
consistent  with  the  stellar  radial  velocity  of  GV  Tau  S  reported  in  Table  4.  The 
stellar  model  used  in  the  fit  is  consistent  with  the  fit  to  the  Na  and  Mg/Al  regions 
(Fig.  2).  The  combined  disk  emission  and  photospheric  absorption  model  {solid 
gray  curve)  is  a  reasonable  fit  to  the  observed  spectrum.  [See  the  electronic  edition 
of  the  Journal  for  a  color  version  of  this  figure .] 


5.1.  Radial  Velocity  Variability  in  GV  Tau  S 

The  stellar  radial  velocity  of  i>lsr  =  9.4  ±1.7  km  s_1  found 
for  GV  Tau  S  differs  from  the  vlsr  observed  in  2001  November 
(—6.2  ±1.5  km  s-1;  Covey  et  al.  2006),  and  the  i;Lsr  in  1999 
December  (3.1  ±  4  km  s  1 ,  see  Table  4;  White  &  Hillenbrand 
2004).  The  vlsr  reported  here  differs  slightly  from  the  6.2  ± 
0.5  km  s_1  velocity  of  the  Taurus  molecular  cloud  in  the  vicinity 
of  GV  Tau,  as  measured  in  the  CO  J  =  1  -0  transition  (Dame  et  al. 
2001)  and  the  dlsr  =  7.0  ±  0.5  km  s_1  measured  for  the  gaseous 
envelope  surrounding  GV  Tau  (Hogerheijde  et  al.  1998). 

A  plausible  explanation  for  the  radial  velocity  variations  ex¬ 
hibited  by  GV  Tau  S  is  that  it  is  a  spectroscopic  binary  with  a 
secondary  that  is  sufficiently  faint  that  it  is  undetectable  in  our 
spectra.  Assuming  that  the  systemic  velocity  is  that  of  the  molecular 
cloud  and  that  the  primary  has  the  stellar  properties  (7j.(r,  M*) 
derived  earlier  (§4.1.1),  the  three  reported  radial  velocities  for  GV 
Tau  S  can  be  accounted  for  if  GV  Tau  S  possesses  a~0. 13  M.  com¬ 
panion  in  a  Porb it  ~  38  day  circular  orbit  (i  =  90,  a  =  0.35  AU).  A 
companion  with  this  mass  at  the  typical  age  of  Taurus  sources 
(~1— 3  Myr),  would  contribute  a  small  fraction  (20%)  of  the 
observed  K  continuum.  Its  A'-band  spectral  features  would  ap¬ 
pear  only  ~2%  deep  in  the  observed  spectrum,  consistent  with 
the  apparent  lack  of  such  features  in  our  observed  spectrum.  The 
observed  radial  velocity  variation  and  the  A'-band  spectrum  are 


Fig.  4. — Observed  is  =  2-0  CO  spectrum  of  GV  Tau  N  {black  histogram) 
also  showing  both  emission  and  absorption  components.  A  model  of  CO  emission 
from  a  rotating  disk  {dashed  line)  is  constructed  from  the  following  parame¬ 
ters:  v  sin  i  (inner  edge)  =  95  km  s”1,  R0/Ri  —  7,  T  —  3000(r/R;)-°-35  K,  E  = 
2.0 {r/Ri)~0  5  g  cm~2,  where  r  is  the  disk  radius.  The  disk  model  is  combined  with 
a  model  of  stellar  photospheric  CO  absorption  {dotted  line)  at  the  same  radial  veloc¬ 
ity  having  the  physical  parameters  Tes  =  4100  K,  log  g  —  4.0,  v  sin  i  —  15  km  s™1, 
and  rK  —  12.0.  At  i>lsr  —  —4.5  km  s-1  (Table  4),  the  combination  {solid gray 
curve)  is  a  reasonable  fit  to  the  observed  CO  feature  in  GV  Tau  N.  [See  the 
electronic  edition  of  the  Journal  for  a  color  version  of  this  figure.] 
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Topocentric  Wavelength  (/im) 

Fig.  5. — High-resolution  (R  —  24, 000)  spectrum  of  GV  Tau  N  (solid  histogram )  with  a  stellar  photospheric  model  fit  (dotted  line )  to  the  absorption  component  that 
appears  to  be  present  within  the  broad  emission  observed  in  the  Na  region  of  the  K  band.  The  Na  absorption  component  is  reasonably  fit  with  the  same  stellar  parameters 
used  in  modeling  the  CO  region  (Fig.  4;  Tes  —  4100  K,  log  g  =  4.0,  vsin  i  —  15  km  s-1,  and  ulsr  =  —4.5  km  s-1)  and  with  veiling  vk  —  15.  The  emission  spectrum 
(dashed  histogram),  obtained  by  subtracting  the  model  absorption  from  the  observed  spectrum,  reveals  broad  emission  in  the  Na  lines  and  possibly  the  Si  line.  [See  the 
electronic  edition  of  the  Journal  for  a  color  version  of  this  figure .] 


also  consistent  with  smaller  companion  masses  and  orbital  sepa¬ 
rations.  A  larger  stellar  mass  for  the  companion  in  a  more  face-on 
orbit  is  ruled  out  since  it  would  be  detectable  in  our  high  signal- 
to-noise  ratio  spectra. 

5.2.  Warm  Molecular  Absorption  in  GV  Tau  N 

Gibb  et  al.  (2007)  previously  reported  absorption  toward  GV 
Tau  S  in  the  HCN  low-,/  R  and  P  lines  near  3.0  /mi,  based  on  ob- 

TABLE  2 

GV  Tau  North:  Radial  Velocity  Measurements 


VLSR 

Order/ Wavelength  Lines(s)  Absorption  Emission  (km  s-1)  Notes 


Order  37,  2.03  gm .  5(2)1 -0  H2  No  Yes  5.9  Line  centroid 

Order  36,  2.11  fim .  Mg  i,  Al  i  No  No 

Order  35,  2.16  gm .  Br7  No  Yes  —17.0  Line  centroid 

Order  34,  2.21  gm .  Na  i,  Si  i,  Sc  i  Yes  Yes  —4.5  ±  4.0  Velocity  fit  to  Na  absorption  cores 

Order  34,  2.22  gm .  5(0)1  -0  H2  No  Yes  8.5  Line  centroid 

Order  33,  2.29  gm .  v  —  2-0  CO  Yes  Yes  —4.5  ±  4.0  Velocity  fit  to  CO  absorption  core 

Order  25,  3.07  gm .  HCN,  C2H2  Yes  No  8.7  ±  0.3  Line  centroid 


servations  made  on  2006  February  17  and  18.  The  HCN  absorp¬ 
tion  was  characterized  by  equivalent  widths  ^0.015  cm-1,  a  radial 
velocity  ^lsr  =  9.2  ±1.9  km  s_1,  and  a  rotational  tempera¬ 
ture  of  T  =  1 15  K.  No  HCN  absorption  was  reported  toward 
GV  Tau  N  in  the  Gibb  et  al.  (2007)  study.  The  wavelength  cov¬ 
erage  of  our  observations  includes  many  of  the  same  lines  that 
were  observed  by  Gibb  et  al.  (2007)  as  well  as  higher  rotational 
lines. 
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Fig.  6. — Resolved  Bry  emission  profiles  for  GV  Tau  S  ( top ),  and  GV  Tau  N 
{bottom).  Zero  velocity  is  defined  as  the  stellar  velocity,  traced  by  the  radial  ve¬ 
locity  of  AT-band  absorption  lines  ( Table  4).  The  vertical  dashed  lines  show  the 
respective  velocity  centroids  of  the  Bry  feature,  while  the  solid  horizontal  lines 
illustrate  the  feature  width  ( FWHM ).  The  velocity  centroid  relative  to  the  stellar 
velocity  is  blueshifted  by  5.7  and  12.5  km  s_I  for  GV  Tau  S  and  N,  respectively 
(Tables  1  and  2). 


At  the  epoch  of  our  observations,  no  HCN  absorption  was 
present  toward  GV  Tau  S.  However,  HCN  absorption  was  de¬ 
tected  toward  GV  Tau  N.  The  absorption  lines  we  detect  toward 
GV  Tau  N  have  equivalent  widths  and  radial  velocities  similar  to 
those  reported  by  Gibb  et  al.  (2007)  toward  GV  Tau  S. 

We  find  a  higher  temperature  (~550  K)  for  the  absorbing  gas 
toward  GV  Tau  N,  but  this  temperature  is  consistent  with  the 
relative  HCN  line  strengths  observed  by  Gibb  et  al.  (2007).  Their 
lower  temperature  was  based  on  analysis  of  only  low-./  lines 
(i.e.,  R1-R6,  P2-P7);  warmer  temperatures  are  indicated  when 
the  higher  J  lines  used  in  our  analysis  are  included.  Consistent 
with  our  results,  HCN  absorption  with  properties  similar  to  those 
reported  here  was  also  detected  toward  GV  Tau  N  (but  not  GV 
Tau  S)  in  mid-infrared  observations  made  with  TEXES  on 
Gemini-North  at  an  intermediate  epoch  (2006  November;  J.  R. 
Najita  et  al.,  in  preparation). 

These  results  may  indicate  that  the  molecular  absorption  to¬ 
ward  both  components  is  variable.  Indeed,  both  GV  Tau  N  and  S 
are  found  to  experience  significant  photometric  variability  in  the 
H  through  M  bands  and  in  the  3.1  /im  ice  feature  ( Leinert  et  al. 
2001).  The  high  temperature  of  the  absorbing  gas  suggests  that 
the  gas  is  located  close  to  the  star  rather  than  in  a  distant  circum- 
binary  envelope  or  the  molecular  cloud.  So  if  variability  accounts 
for  the  difference  in  the  absorption  properties,  we  are  likely  to  be 
observing  the  variations  in  the  near-circumstellar  environment  of 
each  source,  rather  than  the  motion  of  a  distant  absorber  across  our 
line  of  sight.  Such  a  distant  absorber  is  believed  to  account  for  the 
photometric  variations  in  GV  Tau  S  and  perhaps  variations  in  the 
ice  absorption  (Leinert  et  al.  2001). 

Another,  more  mundane,  possibility  is  that  the  components  of 
the  binary  pair  were  misidentified  by  Gibb  et  al.  (2007).  The 
interpretation  that  the  warm  molecular  absorption  is  entirely 
associated  with  GV  Tau  N  would  be  consistent  with  the  system- 
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Fig.  7. —  Spectral  images  of  both  components  of  GV  Tau  showing  extended 
H2  emission.  The  £(0)1  -0  emission  (2.2233  pm,  top)  and  the  £(2)1-0  emission 
(2.0338  pm,  bottom)  are  observed  at  the  positions  of  both  GV  Tau  S  and  GV  Tau  N. 
In  both  images,  a  bright  knot  of  emission  is  also  clearly  present  6.5"  south  of  the 
GV  Tau  S  spectral  trace. 


atically  stronger  3 . 1  //m  ice  feature  measured  toward  this  com¬ 
ponent,  the  strong  silicate  absorption  that  is  detected  only  toward 
GV  Tau  N  (Van  Cleve  et  al.  1994),  and  the  hypothesis  from  the 
literature  that  GV  Tau  N  is  surrounded  by  an  edge-on  disk.  Ad¬ 
ditional  observations  of  the  GV  Tau  system  would  be  useful  in 
resolving  this  issue. 

Setting  aside  for  now  the  possibility  that  the  HCN  absorption 
is  time  variable,  we  can  explore  the  nature  of  the  absorbing  gas 
that  we  observed  based  on  the  stellar  and  circumstellar  properties 
that  we  have  found  for  GV  Tau  N.  As  described  at  the  beginning 
of  this  section,  the  stellar  velocity  is  blueshifted  from  the  en¬ 
velope  velocity  by  ~11  km  s  1 .  This  suggests  that  the  detected 
stellar  component  of  GV  Tau  N  possesses  an  orbiting  companion. 
If  it  does  not,  at  a  velocity  of  1 1  km  s  1  relative  to  the  envelope, 
the  star  would  travel  0.1  pc  in  104  yr,  escaping  its  molecular  en¬ 
velope  on  a  timescale  much  less  than  the  age  of  the  system.  If  GV 
Tau  N  is  then  a  binary,  the  envelope  velocity  is  the  more  appro¬ 
priate  velocity  reference  frame  for  the  system.  Another  possibility 
is  that  the  stellar  light  is  not  seen  directly,  but  rather  in  reflection 
against  material  on  the  far  side  of  the  envelope  that  is  infalling  to¬ 
ward  us  (S.  Strom  2007,  private  communication).  In  either  interpre¬ 
tation,  the  measured  stellar  velocity  is  not  the  systemic  velocity. 


a 


Topocentric  Wavelength  (/xm) 


Fig.  8. — L-band  spectra  of  GV  Tau  N  and  the  telluric  standard,  HR  1412  ( black  histograms )  in  the  topocentric  velocity  frame  where  the  observed  HCN  lines  ( gray 
tickmarks)  are  redshifted  by  35.3  km  s-1.  Spectral  regions  with  poor  telluric  cancellation  have  been  excised  from  the  GV  Tau  N  spectrum.  The  black  tickmarks  denote  the 
expected  positions  of  C2H2  lines  [1/5  and  v-i  +  (^4  +  ^5)+,  top  and  bottom  sets,  respectively).  A  synthetic  model  of  HCN  absorption  in  a  disk  with  temperature  T  =  550  K 
and  column  density  N  —  1.5  x  10l7cm-2  {gray  line )  provides  a  good  fit  to  the  observations.  [See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this  figure.'] 


TABLE  3 

HCN  Lines  in  GV  Tau  N 


Line  ID 

'West 

O m) 

^topo 

Om) 

^LSR 

(km  s-1) 

Percent  Transmission 

Equivalent  Width 

(A) 

pn . 

.  3.0508033 

3.0511622 

8.73 

0.75 

0.12 

pn . 

.  3.0537670 

0.35 

0.21 

P13 . 

.  3.0567549 

0.20 

0.12 

P14 . 

.  3.0597676 

3.0601268 

8.65 

0.98 

0.14 

PI  5 . 

.  3.0628043 

3.0631660 

8.86 

0.70 

0.15 

P16 . 

.  3.0658660 

3.0662214 

8.21 

0.85 

0.19 

P17 . 

.  3.0689521 

0.25 

0.24 

P18 . 

.  3.0720633 

0.45 

0.15 

P19 . 

.  3.0751988 

0.65 

0.20 

P20 . 

Pll,  14,15, 16 . 

.  3.0783594 

3.0787095 

7.56 

8.7  ±  0.3 

0.95 

0.14 
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Fig.  9. — HCN  absorption  lines  in  GV  Tau  N  in  spectral  regions  with  atmo¬ 
spheric  transmission  >70%  (P 1 1 ,  P 14,  P 1 5,  and  P 1 6  transitions,  thin  lines)  aver¬ 
aged  together  ( thick  black  line)  and  fit  with  a  Gaussian  profile  {dashed  line).  The 
velocity  centroid  and  1  o  error  of  the  Gaussian  fit  to  the  average  profile  is  8.7  ± 
0.3  km  s-1  (^lsr)-  [See  the  electronic  edition  of  the  Journal for  a  color  version  of 
this  figure .] 

As  also  described  at  the  beginning  of  this  section,  the  radial 
velocity  of  the  HCN  absorption  is  marginally  different  from  that 
of  the  GV  Tau  molecular  envelope  by  1.7  ±  1.0  km  s  '.  The 
small  velocity  difference  allows  for  the  possibility  that  some  of  the 
absorption  arises  in  gas  that  is  infalling  with  respect  to  the  enve¬ 
lope.  Higher  spectral  resolution  observations  (e.g.,  with  TEXES; 
J.  R.  Najita  et  al.,  in  preparation)  can  better  determine  whether  this 
is  the  case. 

In  the  meantime,  we  can  consider  two  possible  scenarios  for 
the  origin  of  the  HCN  absorption:  (1)  some  fraction  of  the  absorp¬ 
tion  arises  in  gas  infalling  from  the  molecular  envelope  toward  the 
star  at  ^2  km  s_1  and/or  (2)  the  absorption  arises  in  warm  gas  in 
orbit  around  the  star. 


5.2.1.  Origin  in  an  Envelope 

Since  the  bolometric  luminosity  of  the  system  is  7-9  Lq 
(Kenyon  &  Hartmann  1995;  Furlan  et  al.  2007)  and  the  stellar 
contribution  (photospheric  and  stellar  accretion;  §§4.1.2  and  4.2) 
to  the  luminosity  is  ~1.0  Lq,  the  residual  disk  and  envelope  ac¬ 
cretion  isTacc  ~  6-8  Lq.  The  large  residual  accretion  luminosity 
compared  to  the  stellar  accretion  luminosity  for  GV  Tau  N  in¬ 
ferred  from  the  Br7  emission  (~0.2  ZQ;  §  4.2)  implies  that  the 
disk/envelope  accretion  rate  is  much  larger  than  the  stellar  accre¬ 
tion  rate  for  GV  Tau  N,  similar  to  the  results  obtained  for  Class  I 
sources  by  White  &  Hillenbrand  (2004).  For  accretion  in  a  disk 
inward  to  a  radius  of  R\n,  the  disk  accretion  luminosity  Ld  = 
GM^Md/2R[n  of  7  Lq  corresponds  to  a  disk  accretion  rate  of 
Md  ~  1CV6  Mq  yr_1  forR;n  =  2.5  R@  and  a  stellar  mass  of  M*  = 
0 . 8  Mq  .  In  comparison,  the  spectral  energy  distributions  of  most  of 
the  Class  I  sources  in  Taurus  are  fit  with  envelope  accretion  rates 
up  to  an  order  of  magnitude  larger  (Furlan  et  al.  2007).  A  lower 
accretion  rate  ^10~6  Mq  yr~'  is  more  consistent  with  the  weak 
envelope  emission  and  outflow  activity  associated  with  GV  Tau. 

For  the  estimated  mass  of  the  stellar  component(s)  of  GV  Tau  N 
(0.8  Mq),  the  envelope  infall  velocity  tv  =  (2  GMJr) 1/2  is  expected 
to  be  ~2  km  s  1  at  distances  of  360  AU.  In  contrast,  models  of 
gas  in  collapsing  protostellar  envelopes  suggest  that  the  observed 
HCN  temperature  of  ~500  K  can  be  achieved  only  close  to  the 
star,  within  2  AU  given  the  accretion  luminosity  of  the  system 
(Ceccarelli  et  al.  1 996,  their  Fig.  4),  i.e.,  well  within  the  distance  of 
360  AU  that  would  be  inferred  for  the  infalling  gas  based  on  its 
velocity  relative  to  the  molecular  envelope.  It  therefore  appears 
unlikely  that  most  of  the  absorption  could  arise  in  an  infalling  en¬ 
velope  unless  there  are  additional  heating  processes  for  the  in¬ 
falling  gas  beyond  those  considered  in  Ceccarelli  et  al.  (1996).  If 
other  heating  processes,  such  as  oblique  shocks,  can  heat  distant 
material  to  ~500  K,  the  warm  molecular  absorption  that  we  ob¬ 
serve  could  also  arise  in  a  distant,  noninfalling  envelope.  How¬ 
ever,  the  lack  of  molecular  absorption  observed  toward  GV  Tau  S 
argues  against  the  gas  originating  from  a  distant  envelope  which 
would  surround  both  GV  Tau  N  and  S. 

An  alternative  scenario  is  that  GV  Tau  N  is  a  single  star  and  its 
measured  stellar  velocity  is  the  appropriate  systemic  velocity.  In 
this  case,  at  a  distance  2  AU  from  the  star,  the  radial  component 
of  the  infall  velocity  would  be  (GMJr)m  ~  20  km  s_1  (Cassen 
&  Moosman  1981)  for  the  estimated  stellar  mass  of  GV  Tau  N. 
This  is  similar  to  the  13  km  s_1  velocity  difference  between  the 
HCN  absorption  and  GV  Tau  N  stellar  photosphere  ( Table  4).  If 
GV  Tau  N  is  a  single  star,  it  is  then  moving  relative  to  the  mo¬ 
lecular  envelope  and,  as  already  discussed  above,  it  will  travel 


TABLE  4 

Radial  Velocities  for  GV  Tau  and  IRS  46 


System 

Component 

«>LSR 

(km  s-1) 

Reference 

GV  Tau . 

Cloud 

6.2  ±  0.5 

Covey  et  al.  (2006) 

Envelope 

7.0  ±  0.5 

Hogerheijde  et  al.  (1998) 

GV  Tau  S . 

Star 

3.1  ±  3.8 

White  &  Hillenbrand  (2004) 

-6.2  ±  1.5 

Covey  et  al.  (2006) 

9.4  ±  1.7 

This  paper 

GV  Tau  N . 

Star 

-4.5  ±  4.0 
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~0. 1  pc  in  <  104  yr.  It  then  seems  unlikely  that  the  star  and  the 
envelope  are  physically  associated,  which  makes  it  difficult  to 
account  for  the  origin  of  the  infalling  gas  in  this  scenario. 

5.2.2.  Origin  in  a  Circumstellar  Disk 

Another  possibility  is  that  the  HCN  absorption  arises  in  the 
heated  atmosphere  of  a  gaseous  disk  (see  also  Gibb  et  al.  2007). 
Models  of  the  chemistry  of  the  inner  regions  of  circumstellar 
disks  predict  that  HCN  will  be  abundant  at  radial  distances  of  a 
few  AU  (Markwick  et  al.  2002).  Models  of  the  thermal  structure 
of  disk  atmospheres  further  predict  that  the  disk  surface  will  reach 
temperatures  ^  500  K  at  distances  of  a  few  AU  (e.g.,  Glassgold  et  al. 
2004).  If  the  disk  in  the  GV  Tau  N  system  is  viewed  edge-on,  as 
hypothesized  in  the  literature  ( Menard  et  al.  1993;  Van  Cleve  et  al. 
1994),  it  seems  plausible  that  the  warm  HCN  absorption  could 
arise  in  the  disk  atmosphere  seen  in  absorption  against  the  warmer, 
inner  region  of  the  disk  that  produces  the  Z-band  continuum. 

The  rarity  of  such  a  line  of  sight  through  a  disk  may  also  ex¬ 
plain  another  unusual  characteristic  of  GV  Tau:  the  low  CO2  ice 
absorption  optical  depth  toward  the  source  compared  to  its  silicate 
optical  depth.  The  data  of  Furlan  et  al.  (2007)  show  that  the  peak 
10  //m  silicate  optical  depth  is  ~9  times  the  peak  optical  depth  of 
the  15.2  p m  CO2  ice  absorption.  All  other  Taurus  Class  I  sources 
in  the  Furlan  sample  have  rsiiicate/rco2  ~  4,  with  the  exception 
of  DG  Tau  B,  another  possible  edge-on  source.  A  high  ratio  of 
Tsiiicate/Tco2  may  occur  if  a  substantial  fraction  of  the  silicate  op¬ 
tical  depth  is  produced  in  a  warm  disk  region  in  which  CCB  ice 
has  sublimated. 

GV  Tau  N  would  be  a  spectroscopic  binary  in  this  scenario,  in 
order  to  account  for  the  difference  in  radial  velocity  between  the 
stellar  photosphere  and  the  warm  molecular  absorption.  Because 
the  warm  HCN  absorption  is  approximately  at  rest  relative  to  the 
molecular  envelope  velocity,  in  this  scenario  it  would  arise  in  a 
circumbinary  disk  seen  edge-on.  In  contrast,  the  CO  overtone 
emission,  which  is  found  to  share  the  radial  velocity  of  the  stellar 
photosphere,  arises  from  the  circumstellar  disk  of  the  primary 
component. 

5.3.  Comparison  to  IRS  46 

Strong  absorption  by  gaseous  warm  organic  molecules  appears 
to  occur  rarely  among  low-mass  young  stars.  GV  Tau  and  another 
source,  IRS  46,  are  the  only  systems  in  which  such  absorption  has 
been  reported  to  date.  IRS  46  is  a  low-mass  YSO  with  a  spectral 
energy  distribution  that  indicates  either  a  Class  I  source  or  a  Class  II 
source  viewed  nearly  edge-on  (Lahuis  et  al.  2006).  IRS  46  is  similar 
to  GV  Tau  N  in  that  it  also  shows  Bry  and  2.3  /jm  CO  overtone 
emission  (see  YLW  16B  in  Doppmann  et  al.  2005)  and  strong 
absorption  by  gaseous  warm  organic  molecules  in  the  L  band  and 
at  mid-infrared  wavelengths  (CO,  HCN,  C2H2,  CO2;  Lahuis  et  al. 
2006).  The  temperature  (~400  K)  and  equivalent  width  of  the 
HCN  absorption  are  similar  to  that  observed  toward  GV  Tau  N. 
Unlike  GV  Tau  N,  the  warm  HCN  absorption  in  IRS  46  (j'i.sr  = 
—20  km  s~ 1 )  shows  a  large  velocity  offset  from  the  surrounding 
molecular  cloud  (i>lsr  =  4.4  km  s~' ;  Table  4).  In  contrast  to  the 
situation  for  GV  Tau  N,  the  stellar  radial  velocity  of  IRS  46  is  un¬ 
known.  If  the  system  is  at  the  cloud  velocity,  the  warm  molecular 
absorption  is  blueshifted  from  the  system  by  24  km  s-1.  Lahuis 
et  al.  (2006)  speculated  that  the  warm  molecular  absorption  in 
IRS  46  arises  in  the  inner  region  of  a  circumstellar  disk,  possibly 
at  the  footpoint  of  an  outflowing  disk  wind  in  order  to  account 
for  the  presumed  blueshift  of  the  absorption.  This  scenario  is  con¬ 
sistent  with  the  interpretation  of  the  spectral  energy  distribution  as 


indicating  a  disk  system  viewed  edge-on.  Our  “disk  origin”  sce¬ 
nario  for  GV  Tau  N  is  similar  to  the  explanation  given  by  Lahuis 
et  al.  (2006)  for  the  HCN  absorption  in  IRS  46.  In  both  cases,  the 
warm  molecular  absorption  arises  in  a  disk,  although  there  is  no 
outflow  component  to  the  flow  in  the  case  of  GV  Tau  N. 

6.  SUMMARY  AND  CONCLUSIONS 

We  have  used  spatially  resolved  K-  and  L- band  spectroscopy 
to  characterize  the  stellar  and  circumstellar  properties  of  the  GV 
Tau  binary  system.  We  find  that  GV  Tau  S  is  a  radial  velocity  var¬ 
iable,  possibly  the  result  of  an  unseen  low-mass  companion.  The 
radial  velocities  measured  here  and  in  the  literature,  when  com¬ 
bined  with  the  apparent  absence  of  the  spectral  signature  of  a  com¬ 
panion  in  ourK-band  spectra,  are  consistent  with  a  companion  with 
M*  <  0.15  M.  and  a  <  0.35  AU.  Further  spectroscopic  monitor¬ 
ing  of  this  source  would  be  useful  to  confirm  our  interpretation 
and  to  better  constrain  the  companion  properties. 

The  other  component  of  the  binary  system,  GV  Tau  N,  is  found 
to  show  strong  CO  overtone  emission,  strong  A-band  veiling  (ip  = 
12-15),  and  Bry  emission  with  a  blueshifted  emission  centroid; 
signatures  consistent  with  a  young  actively  accreting  star.  From 
the  presence  of  apparent  stellar  absorption  features  in  the  2-0  CO 
overtone  bandhead  and  the  2.2  pm  Na  emission  features,  we  infer 
a  late  spectral  type  (T^g-  ~  4100  K)  and  a  stellar  radial  velocity  of 
r’LSR  =  —4.5  km  s_1.  The  warm  ( T  =  550  K)  HCN  absorption 
found  toward  GV  Tau  N  in  the  Z-band  spectrum  of  this  source  is 
offset  in  velocity  by  ~  1 3  km  s~ 1  relative  to  the  star  and  ~2  km  s~ 1 
relative  to  the  GV  Tau  molecular  envelope.  The  large  velocity  of 
the  star  relative  to  its  molecular  envelope  would  cause  it  to  escape 
the  molecular  envelope  on  a  timescale  much  shorter  than  its  age 
unless  the  star  is  a  spectroscopic  binary. 

The  small  redshift  (1.7  ±1.0  km  s~ 1 )  of  the  warm  molecular 
absorption  relative  to  the  molecular  envelope  velocity  (assumed 
to  be  the  systemic  velocity)  suggests  that  most  of  the  absoiption 
arises  in  a  disk  atmosphere  viewed  close  to  edge  on,  although  an 
origin  in  the  molecular  envelope  is  not  ruled  out  completely.  This 
interpretation  can  be  tested  with  further  observations  of  GV  Tau 
N.  Observations  at  higher  spectral  resolution  can  better  constrain 
the  radial  velocity  of  the  absorbing  gas  relative  to  the  systemic 
velocity.  Observations  that  explore  the  molecular  abundances  of 
the  absorbing  gas  may  also  be  useful  in  testing  a  disk  origin  for 
the  absorption. 
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